PolarPlot[(1 + Sin[9 x]) (1 + Sin[x]) (1 + 0.03 Sin[9*5 x]) (1 + 0.04 Sin[9*33 x]), {x, 0, 2 Pi}]